Towards deterministic equations for Levy walks: the fractional material derivative 
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Levy walks are random processes with an underlying spatiotemporal coupling. This coupling 
penalizes long jumps, and therefore Levy walks give a proper stochastic description for a particle's 
motion with broad jump length distribution. We derive a generalized dynamical formulation for 
Levy walks in which the fractional equivalent of the material derivative occurs. Our approach will 
be useful for the dynamical formulation of Levy walks in an external force field or in phase space for 
which the description in terms of the continuous time random walk or its corresponding generalized 
master equation are less well suited. 
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Anomalous diffusion processes are characterized by 
deviations from the traditional linear time dependence 
(x 2 (t)) = 2Kt of the mean squared displacement in the 
force-free limit. In particular, one distinguishes subdif- 
fusion (0 < a < 1) and superdiffusion (a > 1) for the 
wide class of systems displaying a power-law anomaly 
(x 2 (t)} — 2K a t a /r(l + a); here, K a is a generalized 
diffusion constant [Q, ||. A versatile framework for the 
description of anomalous diffusion are continuous time 
random walks (CTRWs), which define a random walk 
that is governed by two probability density functions 
(pdfs), the jump length and waiting time distributions 
\(x) and ij)(t) from which the jump length x and the 
waiting time t of each jump are drawn |3[|. Although 
the stochastic formulation of the CTRW fully defines the 
random process and leads to the closed integral equation 
for the pdf of the particle's position P(x, t) in terms of 
\(x) and tp(t), its mathematical handling gets awkward 
as soon as non-natural boundary conditions, the pres- 
ence of external force fields, or the description in phase 
space are considered. The same complication holds true 
for the formulation in terms of generalized master equa- 
tions, which are equivalent to CTRWs with uncorrclated 
X(x) and ip(t) Q. In such cases, the corresponding deter- 
ministic equations of the generalized Fokker-Planck type, 
in which the drift terms occur explicitly and which can be 
attacked with the standard mathematical tools, render a 
much more amenable description. To find such equations 
for anomalous transport statistics has been a focal point 
in stochastic systems studies ||. 

For subdiffusion processes, a complete framework of 
transport equations has been established, namely the 
fractional Fokker-Planck and Klein-Kramers equations 
||, 0, fjj . These are natural generalizations of their Brow- 
nian counterparts, and their solution exists whenever 
the solution of the corresponding regular Fokker-Planck 
equation exists, as they correspond to a subordination of 
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the analogous normal stochastic process [|[ 0, |[ ||. 

The description of supcrdiffusive processes within 
the same framework is still far from being completed. 
Whereas Levy flights in the absence of an external force 
display a Levy stable pdf and therefore a diverging mean 
squared displacement (and thus could apply only to 
rather exotic physical processes) jj], |[ |3| [|, [rJ, Levy 
walks (LWs) give a proper dynamical description in the 
superdiffusive domain. The temporal and spatial vari- 
ables of LWs are strongly correlated, their steps being 
governed by a joint distribution ip(x, t), in which waiting 
time and step length pdfs, ip(t) and X(x) are no longer 
independent. In LWs, the occasional long jumps which 
are typical for Levy flights are penalized through the 
introduction of a time cost. This spatiotemporal cou- 
pling can be achieved, in the simplest case, by the choice 
iMx,t) = ^ip(t)S(\x\ — vt), i.e., by a constant velocity 
p|. Such a model arises naturally when describing some 
limiting cases of molecular collisions jflj]; its close rela- 
tives are reasonable candidates for describing turbulent 
dispersion A question of fundamental interest is 

therefore the formulation of LWs in terms of determin- 
istic equations. Whereas previous approaches [|l3| Q in 
terms of fractional Klein-Kramers equations could repro- 
duce lower order moments of an LW, they were hampered 
by the fact that they could not describe the full pdf. The 
main complication on this way is the fact, that the overall 
LW process cannot be immediately considered as subor- 
dinated to a Wiener one (or to a simple random walk); 
however, as we proceed to show, it is exactly the strong 
correlation of the temporal and the spatial aspects of LWs 
which makes it possible to provide a description based on 
a process subordinated to a simple two-state Markovian 
process. In the present work, we derive the exact de- 
terministic evolution equation for LWs which holds both 
for the free motion and in a constant force field. In the 
following, we use rescaled quantities and concentrate on 
the one dimensional case. 

We first define a two-state Markovian random process 
describing the velocity switching and then proceed to 
generalize it to two different domains of LWs. Thus, let 
us denote by P+ and P_ the probabilities to move to the 
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right or to the left, respectively. Probability conservation 
demands that P + + P_ = 1. Moreover, for simplicity we 
assume that the absolute value of the velocity of motion 
to the right and to the left is 1. Within the rate descrip- 
tion, for a symmetric case, the probabilities P± satisfy 
the differential equation 



dt 



P+-P^ 



(1) 



Eq. (|l|) can be readily solved: Taking P_ = 1 — P + we get: 
jj- t P+ = 1 — 2P + . The equilibrium situation corresponds 
to P + = P_ = 1/2, and the relaxation to this equilibrium 
from the initial condition P + = 1 is exponential, 



P± = h ± |exp(-2t). 



(2) 



Let us concentrate first on the switching process de- 
scribed by this equation: It is an alternating random 
process with the waiting time pdf ip(t) = cxp(— t). At 
each "tick" the state is changed from +1 to —1 and back. 
P+(t) is then the probability that at time t, the state of 
the system is "+1", i.e. that the overall number of full 
steps (changes of sign) was even. The Laplace transform 
of this probability is 



P +( u ) = ^2x2n(u) = 



71 = 



u [1 + ip{u)} ' 



(3) 



A similar expression for P_ involves the summation over 
the odd numbers of steps. In Eq. (j|), the X2n(t) denote 
that the walker performs 2n direction changes within an 
overall waiting time t. The latter is denoted by ^(t) = 
1 — Jq ili{t)dt with Laplace transform i^{u) = u~ l —ip(u) ju 
pf. For our exponential function (ip(u) = 1/(1 + u)), the 



result becomes P+ (u) = j- + 



2m ' 2(2+ji) 



, which is exactly the 



Laplace transform of Eq. (||). 

Consider now a a long-tailed waiting time pdf ip(t) 
t~ 1 ~ a of the explicit form |15| 



I 



1 



< a < 1 



(4) 



in Laplace space. This specific form has the following 
origin due to subordination |jq |. In a system whose re- 
laxation in its operational time is given by an exponential 
4>{t) = cxp(— r), we can introduce a coarse-graining in 
which the operational time is divided into intervals At 
(and At taken as a new time unit). Assume that the du- 
ration of the physical time interval corresponding to At 
is given by a one-sided Levy distribution. The duration of 
the physical time corresponding to the interval r is then 
a convolution of n = t/At such distributions, and its 
Laplace transform is exp(— nu a ). Averaging over n we 
get V'( u ) ~ / ex P( — nu a ) exp(— n)dn, exactly reproduc- 
ing Eq. (^). As an example we can explicitly determine 
ip(t) for a = |: ipy 2 (t) = (7rf)~ 1/2 ~ e*erfcVt, with the 
asymptotic behaviors t~ x / 2 (t <C 1) and t~ z / 2 (t 3> !)■ 



With the waiting time pdf (Q), Eq. ([[]) generalizes to 
the fractional form 



(5) 



where qD\ a = 4foD t a , and oD t a is the fractional 
Riemann-Liouville integral operator defined in terms of 



i 



r(a) Jo 



dt' f{t'){t-t') 



f\a-l 



(6) 



with the convenient property / °° e ut oD t a f(t) — 
u~ a f(u) (H 0. In Eq. (||), the fractional derivative on 
the rhs describes a process which is subordinated to the 
simple exponential switching under the operational time 
given by the Levy stable waiting time pdf. To see this let 
us compare two solutions: one using the "CTRW'-time 
and another solving Eq. (||) directly: 

From Eq. (||), with the initial conditions P+(0) = 1 
and P- (0) = 0, we recover upon Laplace transformation 



uP+ - 1 = -M 1 ~ a P + 



uP- = u^P, - u 1_a P_ 



from the second equation, P_ = — Vj Pf, and therefore 
we find by insertion into the first: 



P + 



1 



2u + u° 



and P_ 



1 



2u + u a+1 ' 



(7) 



It is easy to verify that the same result is obtained by 
combining Eq. (||) with Eq. (||). Eq. (0) describes the 
kinetics of moving to the left and to the right. 

We now combine the purely temporal results for P± 
with the drift invoked by a constant velocity, distinguish- 
ing between two different cases. The ensuing propagator 
P of the associated symmetric random walk is combined 
from a superposition of two realizations of the switching 
process, taking place with the rate of 1/2 each, one in 
which the first step goes to the right, and one in which 
it goes left. This introduces an additional factor 1/2 in 
all the following equations. 

(i) Ballistic regime. In the Markovian case, the combi- 
nation of the process (||) with a velocity of magnitude 1 
introduces the material derivatives d+ = M- ± -§-. View- 

111 at ox 

ing now P± as functions of x and t, the evolution equation 
for P±(x,t) result: 



d±P±^\{P^-P±). 
Both together produce the telegrapher's equation 



1 P . 
dt 



^p = *Lp, 

dt 2 dx 2 



(8) 



(9) 



also known as Cattaneo equation [Q, for the quantity 
of interest, the propagator P = P + + P_. The Cattaneo 
equation describes a process which at short times behaves 
ballistically, (x 2 (t)) ~ t 2 , and at long times exhibits nor- 
mal diffusion, (x (t)) ~ t. This derivation was based on 
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the material derivatives d± , whose Fourier-Laplace trans- 
form is u ± ik. We now demonstrate that we reproduce 
exactly the propagator of an LW if we assume (ad hoc) 
that the corresponding fractional material derivative is 
defined through 



T x C t {d^ a f{x, t)} = {u± iky- a f(k, u). 



(10) 



This choice is motivated by the fact that a waiting time 
is still coupled to a walk of length x — t, and that for 
anomalous transport processes the Fourier-Laplace space 
is the natural basis to introduce generalizations. Thus, 
we obtain 



d±P± = U 1 ± - a (r=-r 



(11) 



where the fractional material derivatives are to be inter- 
preted in terms of (|l0l). We solve Eq. under the 
initial condition P + (x,0) = P_(.t,0) = 5(x)/2 so that, 
when introducing the propagator P = P + + P_ and its 
counterpart Q = P + — P_, we have P(x,Q) = 5(x) and 
Q(x,0) — 0. With the abbreviations A+ = u + ik and 
A_ = = u — ik, Eqs. ( pi] ) can be rewritten in terms 
of Fourier-Laplace transformed P and Q, as 



The solution for P reads: 

XIXT 1 + XIXT 1 + XT 1 + A" -1 



P = 



A" + A" + 2\%\° 



(12) 



(13) 



Let us show that this is an exact expression for the 
LW with the waiting time pdf (Q). To this end, note 
that within the CTRW the propagator is obtained as 
P(k,u) = #(k,u)/[l — ip(k,u)] in the case of spatiotem- 
poral coupling with tfj(x, t) — \ [&{ x — t) + $( x + t)] V'W 
and V(x,t) = ±[6(x-t)+6(x + t)]v{t) @. Conse- 
quently, we find 



ip{u, k) = — [ip(u + ik) + ip(u — ik)] 

and an analogous expression for ^(u, k), such that we 
arrive at the Fourier-Laplace form of P: 



P = 



[l-y;(A + )]/A + + [l-y,(A_)]/A. 
2 - V(A+) - V(A_) 



(14) 



With ijj(u) given by Eq. (jfy, Eq. (|13|) is reproduced 
and we have shown that Eq. (|12|) with the defi- 
nition JH]) describes an LW. From the representa- 



tion (114), we find the Laplace space form (x 2 (u)) = 



2 (u a + 1 - a) I (u 3 + 



,3+Q 



j of the second moment, from 



which we obtain the limiting behaviors (x 2 (t)) ~ t 2 for 
t < 1 and (x 2 (t)) ~ (1 - a)t 2 for t > 1, i.e., a mere 
decrease in the amplitude of an overall ballistic process: 
the memory which we introduced by the long-tailed form 
of tp leads to an extreme persistence in a given direction 



on all time scales. There is no turnover to a process with 
a smaller exponent of t as in the Cattaneo case. 

(ii) Subballistic regime. Let us now compare this with 
the better known case of the subballistic domain. We 
again follow our above obtained recipe of formulating 
two equations for the direction-switching process with 
the waiting-time distribution of interest, and then change 
the time derivatives for the material ones. There is an 
heuristic way immediately leading to the equations: in 
the Laplace representation, the system of equations 



uP+-l = \f(u){-P++P_) 
UP- = i/(u)(P+-P_) 



(15) 



leads to the solution P + = 2 / n ^Yiv"\\ ■ Noting that ac- 

T u(u+2/(m)) o 

cording to Eq. (||), this is {u [1 + ^(w)]} -1 , we find the 
relation ip(u) = i-|_ t[ //( M ) ■ H we want a function behav- 
ing for small u as ip(u) ~ 1 — u — it 1+/3 (i.e., one with 
a = 1 + [3 being in the interval between 1 and 2) we have 
to choose f(u) = l + u a , which produces ip(u) — prj^rr^ ■ 



This corresponds to the equations 

;P ± 4(l+o<) (P T -P ± 



d_ 
~dl' 



(16) 



for the alternating process, and to the equations with 
material derivatives 



d±P± = i(l 



4)(P T -P ± 



(17) 



for the Levy walk. Using our formal rules, we get in the 
Fourier-Laplace representation: 



1 -x ± {p±q)- 1 - = t\(i + \ p + )Q 1 



from which we find the result: 



P 



2 + A^ + X't + A + + A_ 
2A+A_ + A+(l + X _) + A_(l + X p + ) ' 



(18) 



This corresponds exactly to Eq. (Jl4|) for the new tp and 
again corroborates the recipe to generalize d± to the frac- 
tional material derivatives d± in the Fourier-Laplace do- 
main. The second moment of this process is obtained 
as (x 2 (u)) = 2{u + flu 13 ) I (u 4 + u 3+l3 + u 3 ), produc- 
ing the limiting behaviors (x 2 (t)) ~ t 2 for t -C 1 and 
(x 2 (t)) ~ 2/3t 2 " /3 /r(3 - 0) for t > 1, i.e., a transition 
from initial ballistic to terminal subballistic superdiffu- 
sive behavior, in analogy to the CTRW result [lSf| . 

Let us now discuss the coordinate-time representation 
of the fractional material derivative. Using the standard 
relation F(u + b) '= e~ bt f(t) of the Laplace transforma- 
tion we obtain after some steps 



d%P(x,t) = Q D?P(x±t,t). 



(19) 



The fractional material derivative, this is, generalizes the 
standard material derivative, d±P{x,t) — 4iP(x±t,t) = 
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(Jj ± gj) P for a = 1, through the introduction of the 
standard Riemann-Liouville operator acting on the entire 
right hand side. 

For subdiffusion, the major advantage of the fractional 
dynamical equation formulation is in the possibility to 
easily generalize to situations with an external force field, 
which led to the fractional Fokker-Planck equation ||, 
^| . Here, we start with incorporating a constant external 
force. To this end, let us consider the physical realization 
of the walk in a splitting flow: in the upper half-plane, 
the particle moves to the right; i>+ = Vq for y > 0, in the 
lower half-plain it moves to the left: v~ — —vq for y < 
0. The motion in the y-direction dictates the waiting- 
time distribution. If it is simple diffusion, the overall 
process is a Levy- walk with a= 1/2 Imagine now, we 
have a force acting in the x-direction. The force causes 
sliding of the particle with respect to the flow, so that 
now w+ = vq + fif for y > and v~ = —vq + [if for 
y < 0. This corresponds to changing d± to the constructs 
corresponding to 

dl ± P(x, t) = D?P(x + \nf ± v]t, t), (20) 

whose Fourier-Laplace transform produces 

df,± -» [u + i(nf±v)k] a . (21) 

That is, in both the force-free and the constant force 



cases, we observe some type of generalized d'Alembert 
principle reflecting the 5-coupling of x and t [^p| . 

Translating the dynamic equations ([h]) and (0) with 
the fractional material derivatives into an equation for 
the propagator P produces a rather complicated expres- 
sion. This can in fact be circumvented by a somewhat 
different definition of the fractional operators, as shown 
in Ref. [|T). However, the latter does not allow for the 
incorporation of a bias and is therefore not suited for our 
purpose. 

In our approach we were guided by the equivalence 
between position x and time t in the LW framework, 
enforced by the (5-coupling tl>(x,t) — hS(\x\ — t)i(>(t), 
which could in fact be rewritten in terms of the jump 
length distribution as — t)X(x), with the appro- 

priate long-tailed form for A. This equivalence gives rise 
to the occurrence of the material derivative, in complete 
analogy to the Brownian Cattaneo case. However, in 
the presence of long-tailed temporal correlations of the 
kind ~ i _1 ~ Q , the fractional variant of the mate- 
rial derivative emerges, with its simple representations 
in both Fourier-Laplace and (x, t) domains. This treat- 
ment is amenable for the case of constant external force. 
Whether there is a similar treatment for arbitrary force 
f(x) is not clear at present. 
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